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ABSTRACT 


Modern  Naval  weapon  and  sensor  systems  are  strongly 
influenced  by  the  marine  environment.  Foremost  among  the 
atmospheric  effects  is  ducting  of  electromagnetic  energy  by 
refractive  layers  in  the  atmosphere.  To  assess  the  effect 
of  ducting  or.  electromagnetic  emissions,  the  Navy  developed 
the  Integrated  Refractive  Effects  Prediction  System  (IREPS) . 
Research  at  Naval  Postgraduate  Scnool  (NPS)  has  led  to 
development  of  a  state-of-the-art  model  which  can  be  used  to 
predict  changes  to  the  refractive  profile  of  the  lower 
atmosphere.  The  model  uses  radroriCr.de  data  and  surface 
meteorological  observations  to  predict  changes  in  refractive 
conditions  ana  low  level  cloud/iog  formation  over  18  to  30 
hour  periods.  The  model  shews  som*  skill  in  forecasting 
duct  regions  wh«n  subsidence  ra*es  can  be  specified  *o 
within  +  /-.0015  m/s.  This  thesis  shows  the  applicability  of 
the  NPS  marine  atmospheric  mixed  layer  model  "to  fleet  tac¬ 
tics.  Atmospheric  refractive  affects  on  specific  emitters 
car.  be  predicted  when  model  predictions  are  used  ir.  cor.jur.c- 
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I.  INTPOD OCTION 

Over  the  centuries  military  leaders  have  learned  no 
appreciate  and  take  advantage  of  the  effects  of  the  environ¬ 
ment,  and  when  they  failed  to  do  so,  have  met  with  disaster. 
As  it  approached  the  islands  of  Japan  in  the  thirteenth  cen¬ 
tury,  Ghengis  Khan's  Mongol  invasion  fleet  was  destroyed  by 
a  typhoon.  From  August  to  October  1599  the  Spanish  Armada 
was  beset  by  storm  after  storm  which  resulted  in  the  sinking 
of  many  ships  and  failure  of  Spain's  attempted  invasion  of 
Britain.  On  D-Day,  June  6,  1944,  the  allied  ir.vastion  force 
took  advantage  of  predicted  good  weather  between  storms  and 
made  *he  amphibious  assault  across  the  beaches  of  Normandy. 
In  December  cf  the  same  year  Admiral  Halsey's  Third  Flee^ 
was  caught  preparing  to  refuel  by  an  undetected  typhoon, 
causing  28  ships  to  be  crippled,  156  airplanes  to  b»  lest, 
and  the  destroyers  S  FENCE,  MONAGHAN,  and  HULL  to  be  sunk 
(Nash  1976). 

Modern  naval  warfare  technology  is  providing  increas¬ 
ingly  capable  and  complex,  but  environmentally  dependent 
wea pons/senscr  systems.  Not  only  the  severe  sea  and  weather 
conditions  must  be  predicted  to  retain  advantage  over  adver¬ 
saries,  but  also  environmental  factors  which  enhance  or 


aggrade  a  wide  range  cf  weapon  and  sensor  systems  which  uti¬ 
lize  electromagnetic  (EM)  propagation.  Electromagnetic  fre¬ 
quencies  above  the  HP  (3-30Mhz)  band  can  be  greatly  affected 
bv  atmospheric  refraction. 

A  condition  known  as  ducting  occurs  when  a  refractive 
layers  cause  EM  energy  to  bend  toward  the  ear-1-  h  at  a  rate 
greater  than  or  equal  to  the  earth's  curvature.  Ducting 
occurs  with  certain  critical  vertical  gradients  of  tempera¬ 
ture  and  humidity,  and  can  cause  both  increased  radar  and 
radio  ranges,  and  holes  (gaps)  in  normal  coverage. 

Tactical  advantages  exist  by  knowing  duct  locations. 
T-.sse  advantages  include  being  able  tc  make  realistic  esti¬ 
mates  cf  ESM  detection  and  counter  detection  ranges.  This 
would  enable  commanders  to  make  decisions  concerning  emis¬ 
sion  control  (EMCON)  ,  the  positioning  of  both  air  and  sur¬ 
face  surveillance  assets,  the  altitude  and  flight  profiles 
for  strike  aircraft  to  minimize  detection,  the  placement  of 
electronic  gammers  for  maximum  effect,  and  numerous  other 
tactical  considerations. 


Ducting  commonly  occurs  with  inversions  which  act  as 


aloft  and  cooler,  more  moist  air  below  and  typically  exist 
in  marine  surface  high  pressure  regions.  Extensive  areas  of 
low  level  stratus  clouds  often  delineate  areas  of  duct 
cocurren ce . 

Ductir.g  is  expected  to  be  minimal  near  fronts  and  areas 
cf  convective  cloudiness.  Fronts,  with  their  associated 
upward  motion,  often  dissipate  the  inversions  as  the  whole 
air  column  becomes  mixed.  Areas  of  convective  activity, 
detectible  by  the  presence  of  cumulus  clouds,  are  also  nor¬ 
mally  inversion  free.  Thus  satellite  infrared  (IP)  and  vis¬ 
ible  photography  should  provide  means  to  estimate  ducting 
regi  mes. 

The  Navy  employs  a  microcomputer  based  system,  IHEPS 
(In teqra red  Pefractiv*  Effects  Prediction  System)  (Hitnev, 
1?“9),  to  identify  ducting  conditions  and  :o  assess  the 
effects  or.  various  fleet  EM  emitters.  IREP5  requires  radio¬ 
sonde  data  as  input. 

The  disadvantage  cf  IPEPS  is  -hat  it  predicts  the  duct¬ 
ing  conditions  only  at  the  radiosonde  launch  site  and  only 
at  the  sounding  time.  As  such,  it  is  not  an  1REPS  weakness 
but  more  an  inherent  weakness  in  single  station  assessments. 
The  variabilitv  of  the  atmosphere  makes  projections  in  space 
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an  d  predict  ions  in  tire  of  the  ducting  environment  important 
for  the  tactician.  While  satellite  data  can  provide  some  of 
the  data  necessary  to  make  projections  in  spac<=,  something 
else  is  needed  to  predict  changes  with  time.  Considering 
EMCON,  the  ability  tc  forecast  ducting  for  a  24  hour  period 
may  be  critical.  Under  strict  EMCON  conditions  radiosondes 
cannot  be  launched  because  the  signal  they  emit  while  trans¬ 
mitting  data  could  act  like  a  beacon  to  hostile  forces. 

A  predictive  model  exists  for  changes  in  the  marine 
atmospheric  mixed-layer.  The  micro  computer  based  model  was 
developed  by  the  Environmental  Physics  Sroup  at  the  Naval 
Postgraduate  School  (Davidson  et  al,  1933).  The  model  is 
initialized  with  IREPS  sounding  data,  surface  observations 
of  sea  surface  temperature,  wind  speed,  and  subsidence  at 
the  'inversion. 

0*iiizir.a  :he  HP-9845  (which  all  carriers  currently  have 
cr.  board  for  use  as  part  of  the  IREPS  system)  and  data  in 
IREPS  format,  Brower  (1982)  developed  a  program  to  incorpo¬ 
rate  the  KPS  mixed  layer  model  to  predict  changes  in  Ahs 
refractivity  profile  for  a  33  hour  period.  The  model 
requires  as  inputs:  surface  wind,  sea  surface  temperature, 
and  an  estimated  subsidence  value  as  well  as  current  IREPS 
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sounding  data  to  oredic^  the  mixed  layer  evolution  an 
resulting  changes  in  both  surface  based  and  elevated  due 
heights  up  tc  1230  m  altitude. 

The  purpose  of  t  h  is  thesis  is  tc  evaluate  this  state  c 
-he  art  atmospheric  mixed  layer  model  and  its  application  t 
fleet  tactics.  The  data  utilized  ir.  this  study  ware  fro 
radiosondes  taker,  by  U3S  NISI!"  and  OSS  CORAL  SEA  in  th 
Arabian  Sea  in  the  vicinity  of  2 iV  ,  a  data-poor  area 
during  February  1980. 


40°  E 
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80°  E 


Figure  1.  Regional  Chart  Showing  Sounding  Data  Locations. 
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II.  THE  HOD EL  DESCRIPTION 


A.  MODEL  OVERVIEW 

The  NPS  mixed  layer  model  is  ar.  integrated,  -wo  layer, 
zero-order  model  which  predicts  changes  in  the  marine  atmos¬ 
pheric  boundary  layer  (MASLJ  below  the  inversion.  The  two 
layers  consist  of  the  lower,  well-mixed,  turbulent  boundary 
layer  topped  by  the  relatively  non-turbuier.r  free  atmos¬ 
phere.  The  two  layers  are  separated  by  ar.  inversion  or 
transition  zone  which  is  assumed  to  have  zero  thickness. 
The  predicted  properties  are: 

(1)  height  of  inversion 

(2)  values  of  well-mixed  properties 

(3)  values  of  lumps  at  the  inversion 

(4)  formation  of  clouds /fog  within  the  mixed  layer 
Procedural  uses  are  shown  in  Figure  2,  and  the  equations 
will  be  presented  in  Section  3.  Inputs  to  ‘hs  model  are: 

(1)  Radiosonde  data 

(a)  Vertical  distribution  of  temperature 

(b)  Vertical  distribution  of  moisture 

(2)  Sea  surface  temperature 

(3)  Wind  speed 


(4)  Subsidence 

(5)  Latitude 

(6)  Julian  date 

(7)  Local  tine  of  sounding 

From  the  predicted  properties,  a  refractdvity  profile  is 
produced  using  *he  refractivity  equation  (Section  B) . 
Observed  and  predicted  sea  surface  te^oera-.ure  and  wind 
speed  are  input  parameters.  Subsidence  tares  are  estimated 
from  previous  observations.  Latitude,  Julian  c.ate,  and 


local  times  are  used  to  estimate  incident  solar  radiation. 
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EASIC  MODEL  EQUATIONS 


Once  initial  inputs  have  been  made  ar.d  the  vertical 
s-ructure  digitized  by  fitting  straight  lines  *.o  approximate 
the  vertical  profiles  of  temperature  and  humidify,  Equations 
1 A ,  IB,  and  1C  are  used  to  calculate  the  flux  scaling  param¬ 
eters.  The  calculations  are  based  on  the  air-sea  erri  na¬ 
ture  and  specific  humidity  differences  and  the  wine  speed. 
The  scaling  parameters  are  then  used  to  estimate  surface 
fluxes  of  momentum,  heat,  and  moisture  in  Equations  2A,  23, 
and  2C. 

Equation  3 A  is  used  to  predict  the  change  in  mixed- layer 
properties  by  the  difference  between  surface  and  itvrriior. 
fluxes  scaled  by  inversion  height.  Equation  3B  is  use:  *ro 

predict  changes  in  the  inversion  height  by  adding  the  tlux 

• 

at  the  inversion  scaled  by  the  inversion  jump,  a  measure  of 
the  strength  of  the  inversion,  to  the  subsidence.  Equation 
3C  models  the  change  in  the  inversion  jump  by  combining  the 
effects  of  surface  and  inversion  fluxes  with  subsidence  and 
-he  change  cf  -he  inversion  height  determined  ir.  aqua*  ion 
3B.  This  requires  use  of  the  gradient  of  the  proper- y  above 
the  inversion. 
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Flux  Scaling  Parameters: 


0*  *  Cd'^O  (U) 

T*  =  (0  -  9o)  Ce'^  (13) 

q*  =  Ce^fq  -  qc)  (1C) 

Surface  Fluxes: 

-  (o7?*)  =  o*2  ( 2 a. ) 

-  (W^F*  )  =  U*T*  ( 2B) 

-(WTq*)  =  U*q*  (2C) 

Where: 


C d , Ce  =  stability  dependant  drag  coefficients 
0  =  wind  speed 

U*,T*,q*  =  scaling  parameters  of  vertical  turbulent 
momentum,  temp,  and  moisture  -rar.sfer 

-  (u  '  W ' )  =  downward  turbulent  transfer  of  momentum 

-  (W ' T* )  =  upward  turbulent  transfer  of  heat 

-  (W  'q ')  =  upward  turbulent  transfer  cf  moisture 

Predictive  Equations  (general  form)  : 

dXm/dt  =  (WHc’c  -  W^X *  i  *  f  (R)  )  /  h  (3A) 

-wT'i  =  DX  (dh/dt  -  Ws)  ♦  f  ( R)  ( 3 B ) 

d(DX)/dt  =  Yx  ( dh/dt -Ws)  -  (W»"X'  o-WMC'  i+f  (R)  )  /h  (3C) 


Where: 


Xm  =  ary  well- mixed  property:  9,  9v,  or  q 
9  =  potential  temperature 
q  =  specific  humidity 

W '  X '  =  turbulent  vertical  transfer  (flux)  of  X 
o  subscript  =  "at  the  surface" 
i  subscript  =  "at  the  inversion" 
b  =  mixed-layer  depth  (or  inversion  height) 

DX  =  change  in  X  across  the  inversion:  (Xabove 
Ws  =  subsidence  at  top  of  inversion 
Yx  =  vertical  gradient  of  X  above  inversion 
f  (R)  =  radiation  factor 

Pefractivity  Equation: 

a  =  77.6  P/T  ♦  6  X  10*  ( q  P/T* )  ♦  .157  Z 

Where: 

M  =  modified  r efr activir y 


-  Xm) 


P  =  atmospheric  pressure  in  millibars 


HEIGHT  (m) 


figure  3.  Vertical  Structure  of  9,  q  and  Corresponding  M 
Profile  for  a  Typical  Inversion  (Model 
Simplification  rs  Dashed)  . 

Fig.  3  illustra  *es  the  relationship  between  vertical 
distribution  of  temperature  and  humidity  and  the  modified 
refr activity,  M.  M  profiles  are  used  as  model  output 
because  of  the  ease  with  which  ducting  information  car.  be 
extracted  from  them.  The  top  of  a  duct  corresponds  to  the 
heiqht  above  the  surface  where  '■he  M  value  is  a  minimum. 
The  duct  bass  corresponds  to  the  height  at  which  a  vertical 
line  drawn  downward  from  the  point  of  minimum  M  value  first 
intersects  a  point  of  equal  M  units  or  the  surface. 
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SATELLITE  DATA  IN  EOT  TO  MODEL 


I 

s 


Provision  is  made  in  the  mod«l  for  incorporating  fore¬ 
cast  wind  and  sea  surface  temperature  rhanges  during  the 
prediction  period.  Although  the  mixed  layer  is  sensitive  to 
SST ,  it  is  believed  that  accuracies  of  «■/-  1°C  in  SST  are 
sufficient  for  reasonable  model  accuracy.  Stewart  (198  1) 
reports  the  following  satellite  measurement  capability  has 
been  achieved:  NOAA-6  (VHSR)  infrared  radiometer:  SST  to 

•6°C  (with  no  clouds)  ;  SeaSat  and  Nimbus-7  microwave  radiom¬ 
eter  (SMMP.)  :  SST  to  1.0°C  (with  no  rain,  no  P.FI  (radio  fre¬ 
quency  interference)  ,  and  >600  km  from  land)  . 

Satellite  SST  could  be  used  in  the  flux  scaling  equation 
(13)  to  determine  Qo.  Three  methods  are  currently  used 
aboard  ship  for  SST  determination:  Bathythermographs,  sea 

water  injection  temperature,  and  the  ' bucket /thermometer ' 
method.  None  of  these  yield  the  actual  sea  surface  "skin" 
temperature  which  is  the  relevant  quantity  in  the  determina¬ 
tion  of  stability  and  turbulent  heat  and  moisture  fluxes. 
Therefore,  satellite  measurements  would  not  only  provide 
cloud  descriptions  over  broad  areas,  but  would  provide  a 
means  to  obtain  the  "skin"  temperature. 
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III.  ARABIAN  SEA  CLIMATOLOGY  AND  SYNOPTIC 
- BESCffTPTTUN - 


A.  CLIMATOLOGY 

The  climatology  of  the  northern  Arabian  Sea  is  dominated 
tv  the  summer  (southwest)  and  winter  (northeast)  monsoons. 
The  data  used  in  this  study  were  from  February  during  the 
latter  part  cf  the  northeast  monsoon  which  is  characterized 
by  light  offshore  winds  averaging  5-ID  tr.ots,  and  normal 
Hadley  circulation  with  the  subtropical  iet  located  near  30 
N  at  the  200  mb  level.  Fig.  4  shows  the  average  per  cent  cf 
•rime  that  elevated  ducts  occur  during  February  to  March  in 
tais  region. 


Fiqure  4.  Frequency  cf  Elevated  Duct.  Occurrence  Durina 
Spring,  after  Lamtners  et  al  (1990). 


3.  SYNOPTIC  CONDITIONS 

FNOC  (Fleet  Numerical  Oceanographic  Center)  products 
were  used  to  describe  the  synoptic  situation  during  the  da~a 
period.  They  were:  Surface  Analyses,  500  mb  analyses  with 
Temperature  fields,  300  mb  analyses,  25  0  mb  wind  fields. 
Also,  NOAA-6  satellite  photos  were  used  *■  o  determine  cloud 
coverage  of  the  area  as  well  as  to  confirm  jet  positions. 

The  synoptic  situation  in  the  region  from  6-23  February 
1980  was  characterized  by  a  stable  high  pressure  surface 
ridge  over  the  Arabian  peninsula  and  Arabian  Sea,  with  low 
pressure  troughs  over  the  Fed  Sea  to  the  west  and  Indian 


subcontinent  to 


The  iosi  significant  synoptic 


the  east. 

change  occurred  from  11  to  14  February  as  an  853  mb  low 
moved  across  the  northern  part,  of  the  region.  This  low  was 
northwest  of  the  Persian  Gulf  a  t  30  degrees  North  latitude 
on  11  February  and  tracked  eastward,  arriving  over  northern 
India  on  14  February.  It  was  during  this  brief  period  of 
11-14  February  that  the  minimum  shipboard  surface  pressures 
for  the  entire  period  were  recorded.  The  minimum  recorded 
surface  pressure  was  1012.0  mb  on  14  February. 

Table  I  lists  cloud  coverage  over  the  immediate  area  of 
the  da*a.  Low  level  stratus  clouds  were  confirmed  only  on 
13  February  and  from  19  to  21  February. 


r; 


TABLE  I 

Data  Area  Cloud  Coverage  8-21  February  1980. 


DATE 

CL  COD  TYPE  PRESENT 

h- 

8-11 

clear  with  scattered 

cirrus 

1  1-12 

no  satellite  coverage 

1  3 

low  level  stratus 

14-18 

clear  with  scattered 

cirrus 

19-20 

cirrus  and  stratus 

20-21 

cirrus,  stratus  with 

m  u  1 1  i  1 

« 

*  det  erm 

ined  from  NOAA  5  IR  imag 

er  y 

I  - 

r ' 


Ar.  attempt  was  made  to  derive  subsidence  values  from 
synoptic  scale  wind  fields  by  calculating  divergence  at  var¬ 
ious  levels.  It  was  noted  that  at  upper  levels,  wind  maxima 
or  jet  streams  had  a  dramatic  influence  or  divergence  pro¬ 
files.  Although  the  magnitudes  of  vertical  motion  obtained 
in  this  manner  were  inconsistent  and  highly  variable,  an 
interesting  and  possibly  significant  observation  was  made. 

The  subtropical  jet  associated  with  the  wind  maxima  dis¬ 
cussed  above  is  usually  found  at  the  northern  limb  of  the 
Hadley  Cell  (Fig.  5).  It  is  at  this  northern  limb  where  the 
greatest  downward  motion  or  subsidence  is  found  (Palmen  and 
Newton,  1  965).  Although  this  is  a  feature  of  the  large 

scale  circulation,  it  seems  evident  *hat  the  jet  does  have 
>.n  effect  which  should  be  considered  in  single  s-.ation 
assessments  whenever  possible. 

Changes  ir.  the  subsidence  rates  calculated  by  hindcast¬ 
ing  and  shown  in  Table  II  correlated  closely  with  changes  in 
the  position  and  strength  of  the  upoer  level  jet  stream. 
During  the  period  6-9  February  there  was  a  sready  decrease 
in  subsidence  races  from  -.0070  ro  -.0055  m/s.  During  the 
same  period  the  tail  of  the  25  0  mb  jet  max  moved  steadily 
eastward  from  a  point  due  north  of  the  data  area.  From  6-11 


25 


February  the  300  mb  jet  max  shifted  northward  away  from  the 


oat  a  are  <?  • 


TABLE  II 

Calculated  Subsidence  Fates 


RUN  123455789 

DATE  5-7  6-7  7-8  8-9  8-9  1  1-12  14-  15  21-22  22-23 

METHOD 

H  -  -  30  -.70  -.  65  -.60  -.55  -.55  -.75  -.  31  -.  50 

Q  +  -1.72  ♦  -.39  -.06  -1.17  -.55  -.41  -.  33 

A  -.54  -.46  -1.10  ♦  ♦  -.37  -.33  ♦  -.50 

*  subsidence  rates  ir.  gm/s. 

**  H  hindcast,  Q  specific  humidity,  A  adiabatic 


From  9-12  February,  subsidence  rates  remained  at  -.0055 
iti/s.  A  .so,  r.o  change  in  strength  or  position  of  the  250  mb 
jet  occurred  between  11  and  12  February. 

A  subsidence  rate  of  -.0075  m/s  was  determined  for  the 
period  between  14  and  15  February  which  was  an  increase 
since  the  12th.  The  300  mb  jet  strengthened  and  shifted  to 
the  south  toward  the  data  area  from  12-15  February,  and  the 
250  mb  jet  increased  in  strength  from  14-15  February. 

Curing  the  21-2  3  February  period,  subsidence  rates 
increased  fret  -.0031  to  -.0050  m/s.  During  this  oeried. 


THE  MODEL  PERFORMANCE 


A.  EUCT  PREDICTION 

Each  soundin'?  was  run  with  I  REPS  to  establish  a  baseline 


*o  assess  the  model’s  ability  to  predict  changes  ir.  ducting 
conditions.  It  was  then  possible  to  compare  the  model  pre¬ 
dicted  duct  heights  with  those  observed.  Seven  different 
pairs  of  soundings  were  examined,  with  each  pair  separated 
by  24  hours.  Two  of  the  24  hour  periods  had  2  simultaneous, 
spatially  separated  soundings  at  the  end  of  the  periods.  On 
one  occaisicr.,  one  cf  the  second  soundings  was  ■‘taken  at  a 
location  4  25  nautical  miles  to  the  southeast  of  the  first 
sounding . 

The  observed  winds  at  the  start  and  end  of  the  period 
were  used  and  interpolated  linearly  within  the  period  to 
specify  the  wind  for  the  prediction  period.  Sea  surface 
temperature  at  the  latitude  and  lonoitude  of  the  initial 
sounding  and  at  the  latitude  and  longitude  of  the  verifying 
sounding  were  also  interpolated  within  tr.e  period. 

Subsidence  velocities  were  determined  by  hindcasting 


such  that  the  predicted  inversion  height  agreed  with  that 


observed  in  -he  verifying  sounding.  Predicted  ducts  using 
this  "hi recast "  subsidence  are  compared  with  IREFS  observed 
ducts  in  Table  III  .  The  SMS  errors  of  the  model  are 
approximately  one  half  the  errors  of  persistence.  Persis¬ 
tence  is  considered  to  be  a  "prediction"  of  no  change.  It 
is  frequently  used  as  a  baseline  against  which  predictions 
are  compared  to  determine  if  the  predictive  method  has 
merit.  These  results  indicate  that  when  subsidence  car.  be 
accurately  estimated,  the  model's  performance  is  clearly 
superior  to  persistence  in  predicting  duct  bases  and  tops. 


TABLE  III 

Comparison  of  Observed  Duct  Heiahts  With  Model  Predicted 
Heights  Using  Hindcast  Subsidence  Values. 


RUN 

INITIAL 

VERIF  YIN3 

MODEL 

NUMB 

??.  SOUNDING 

SOUNDING 

PREDICTION 

1 

788-983 

470-  877 

345-740 

2 

788-983 

633- 988 

500-960 

3 

470-877 

46- 504 

0-370 

4 

47-504 

0-135 

0-170 

c 

47-504 

0-  304 

0-158 

6 

to  duct 

5-584 

115-375 

7 

^  2  3- 103  1 

248- 569 

0-495 

8 

arc  0-112 
206-  50  3 

no  i uct 

50-420 

9 

r.o  duct 

0-  258 

0-250 

RMS 

ERRORS: 

PREDICTION 

PERSISTENCE 

TOP 

173 

364 

BASE 

1  11 

214 

THICKNESS 

185 

285 

*  heights  in  meters. 
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Figure  6.  Hcdel  Predicted  Ducts  (dashed)  vs.  Observed  Duc*s 
(solid)  for  the  Period  6-8  February  1980. 

For  a  true  predictive  application,  at  least  3  sequential 
soundings  are  required.  The  first  and  second  sounding 
enable  the  estimation  of  subsidencs.  The  model  is  then  ini¬ 
tiated  from  the  second  sounding  using  the  subsidence 
obtained  from  the  first  and  second.  The  third  sounding  is 
used  to  verify  the  prediction.  Gleason  (1982)  found  that  a 
method  which  is  based  on  the  well  mixed  specific  humidity  is 
an  accurate  method  to  calculate  subsidence.  However,  assum¬ 
ing  persistence  in  subsidence  rates  appears  to  be  the  best. 
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A  steady  decrease  in  subsidence  rates  (Table  II)  during 
the  first  week  of  the  data  period  coincided  with  a  steady 
decrease  in  surface  pressure.  The  upper  level  jet  stream 
shifted  away  from  the  data  area  and  weakened  during  this 
period.  This  suggests  that  judicious  modification  of  per¬ 
sistence-obtained  subsidence  may  be  appropriate  in  some 
cases,  especially  when  a  trend  has  been  observed. 

From  three  periods  having  3  sequential  soundirqs.  Tables 
IV  and  V  shew  the  results  of  predictions.  In  the  first 
period  (Hun  number  3)  ,  an  elevated  layer  became  nearly  sur¬ 
face  based.  The  model  prediction  was  for  a  surface  based 
duct.  In  the  second  period  (Hun  number  U)  ,  a  low  elevated 
duct  became  surface  based.  The  model  prediction  matched  the 
observation  nearly  perfectly.  Ir.  the  third  case  (Hun  number 
9)  ,  a  surface  based  duct  formed  when  there  was  no  duct  ini¬ 
tially.  In  this  case,  the  model  predicted  the  formation  of 
an  elevated  duct.  While  the  predicted  duct  in  the  third 
case  differed  substantially  in  both  heioht  ai.d  thickness 
from  the  observed  duct,  the  model  prediction  was  significant 
because  the  occurance  of  a  duct  was  predicted  from  an  ini¬ 


tial  non-ducting  condition 


For  this  particular  comparison. 


I 

* 

i 

* 

'*• 

*  -• 

i 


the  absolute  errors  are 
more  neaningful  than  rms  errors.  Two  of  three  model  pre¬ 
dictions  are  essentially  on  the  mark  while  all  of  the  error 
is  concentrated  in  the  third.  In  the  persistence  cases,  the 
errors  are  fairly  evenly  distributed. 


TABL3  IV 

Comparison  cf  Observed  Ducts  with  Modal  Predicted  Ducts 
Assuming  Persistent  Subsidenoe. 


PUN 

INITIAL 

VERIFYING 

MODEL 

NUMB  ER 

SOUNDING 

SOUNDING 

PREDICTION 

3 

4  70-  87  7 

46- 504 

0-350 

4 

46-504 

0-  135 

0-140 

9 

no  duct 

0-  258 

490-820 

*  heights  in  meters. 


Vi 


I; 

!  • 


32 


TABLE  7 


Error 

Analysis  of 

Model 

Prediction  vs. 

Persist 

e  n  c  e . 

ABSOLUTE  ERROR 

3  RUN 

RMS  ERROR 

RUN 

MODEL 

PERSISTENCE 

MODEL 

PERSISTENCE 

3 

TOP 

154 

373 

336 

337 

BASE 

46 

424 

284 

246 

THICKNESS 

108 

51 

75 

240 

u 

TOP 

5 

369 

3ASE 

0 

47 

THICKNESS 

5 

3  23 

S 

TOP 

562 

258 

BASE 

490 

0 

THICKNESS 

72 

253 

*  heights  in  meters. 


B.  CLOUD  PREDICTION 

The  cloud  formation  is  predicted  on  the  basis  of  concur¬ 
rent  prediction  of  t'n«=  lifting  condensation  level  (LCL)  and 
the  inversion  height.  If  the  LCL  is  below  the  inversion, 
clouds  occur;  if  it  is  above,  no  clouds  occur.  when  model 
predictions  showed  no  low  level  clouds  formina,  satellite 
imagery  confirmed  there  were  no  low  level  clouds  and  that 
the  sounding  site  either  was  clear  or  showed  scattered  high 
cirrus  clouds.  The  formation  of  low  level  stratus  clouds 
was  predicted  on  3  occasions  over  the  entire  period.  Stra¬ 
tus  clouds  were  predicted  to  form  at  2100  local  time  11  Feb¬ 
ruary.  The  first  satellite  imagery  available  following  this 
time  was  at  0800  13  February,  and  showed  the  presence  of  lew 
level  stratus  clouds  (Fig.  8)  . 
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Fioure  7.  Model  Prediction  Showing  LCL  Droppina  Beneath  the 
Inversion  at  2100  11  Feb. 

Clouds  were  predicted  to  form  at  0500  15  February  (25 
hours  into  the  30  hour  prediction).  In  this  case  the  LCL 
dropped  slightly  below  thc  inversion,  but  less  sha'rply  than 
in  the  other  cases.  Satellite  imagery  from  15  February 
shows  no  clouds  (Fig.  10)  . 

Satellite  imagery  shows  the  presence  of  low  level  stra¬ 
tus  clouds  on  the  2000  pass  19  February  and  increased 
amounts  on  the  0800  20  February  pass.  Although  no  soundings 
wsr<=  available  to  initialize  the  model  prior  to  this  initial 
formation  of  clouds,  on  the  next  available  model  run  clouds 
were  predicted  throughout  the  period  0700  21  February  to 


.  Satellite  imagery  confirm  the  presence  of 


0700  22  February 


Figure  8.  N0AA.6  Satellite  13  Imagery  from  0800  13  February 
Showing  tne  Preset  cs  of  Stratus  clouds. 

both  cirrus  and  stratus  clouds  or  both  the  0800  ar.d  2000 
passes  on  21  February  (Fig.  12?  . 

Although  the  satellite  imager/  analysis  was  by  neces¬ 
sity,  subjective,  some  objertivetv  was  retained  by  confirm¬ 
ing  the  initial  analysis  wi-h  a  second  analysis  by  an 
independent  source.  Both  analyses  aareed  or.  cloud  types. 
Since  the  model  could  be  shown  to  be  incorrect  only  on  one 
occasion,  the  model’s  overall  performance  in  predicting  the 
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Figure  12.  FOAA  6  Satellite  IR  Imagery  3300  21  February. 


Figure  11.  F.cdel  Predic 


400 

200  j- 

! _ i 

0700  1300 


for marion.  at  low  level  clouds/fog  was  sight  out  of  tin® 
cases  correct  or  89%  correct. 

C.  EXAMINATION  CF  MOTEL  SENSITIVITY 

Surface  observation  of  sea  surface  temperature  (SST)  and 
wind  can  be  obtained  with  an  accuracy  of  +/-  1  °C  and  +/-  2 
knots  respectively.  To  test  the  model  sensitivity  to  these 
possible  measurement-  errors,  successive  model  runs  ware  made 
with  SST  1  °C  higher  and  lower  char,  observed.  The  effect  of 
wind  measurement  error  was  evaluated  using  wind  values  2 
Vr.ots  higher  'har.  observed.  The  resulting  duct  predictions 
ware  than  compared  with  ducts  predicted  by  actual  observed 
winds  and  SSTs.  Table  VI  contains  the  tabulated  results. 

Noting  that  IMS  error  in  duct  specifications  due  to  pc-en- 
-ial  351  measurement  error  are  approximately  double  the  RMS 
error  due  to  potential  wind  measurement  error,  one  might 
reasonably  conclude  SST  determination  is  the  more  critical 
of  the  two  parameters.  It  is  imperative  to  recognize  that 
the  wind  has  a  much  wider  range  of  variability  than  the  SST 
and  accurate  specification  of  both  parameters  are  important 
for  good  model  performance. 
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?  A  B  LE  V 


Sensitivity  of  Model  Duct  Predictions  to  SST  and  wind 


Varna t ro  ns. 

RUN 

VERIFYING 

SST 

SST 

WIND 

WIND 

NUMBER 

DUCTS 

♦  1°C 

-1°C 

♦  2  KTS 

-2  KTS 

1 

345-740 

400-300 

3 15-715 

350-310 

330-725 

2 

500-960 

730-  1  ’90 

420-870 

625-1080 

545-975 

3 

0-370 

0-395 

0-3  60 

0-400 

0-365 

4 

0-170 

0-270 

0-110 

0-190 

0-140 

5 

0-  153 

0-210 

0-1  15 

0-170 

0-130 

6 

1  1  5-375 

170-570 

0-80 

25-3  1  C 

0-240 

7 

0-495 

290-430 

3-4  4  0 

0-510 

0-450 

8 

5  0-420 

0-360 

0-3  10 

70-450 

0-370 

9 

0-  250 

0-24  0 

0-290 

0-240 

0-19) 

Combined 

RMS 

error  : 

SST 

DU  CT 

WIND 

114 

top 

5  5 

104 

base 

50 

107 

thickness 

26 

*  he 

ig  V  - 

s  in  mete 

rs. 

The 

model’s  demonstrate 

d  sensitivi 

ty  to  sea  surface 

temperat 

ur  o 

points  to 

■5  h  ®  r.  e 

c-d  for  accurate  SST 

det  ermir.a- 

tior.s. 

cur 

rent  methods  to 

measure  sea 

surface 

temperatur® 

ir.cl  ude; 

mcr.it  or  ir.  g 

sea  war 

?r  injection  tempera 

ture  (20-40 

ft.  beneath  the  surface) ;  utilizing  3T  drops  (get  tempera¬ 
ture  at  atou-  1  ft.);  and  the  bucket  and  thermometer  method 
(get  temperature  of  ton  few  inches).  None  of  these  methods 
yield  the  se?  surface  skin  temperature  which  is  the  quantify 
necessary  to  determine  the  surface  fluxes.  Thes®  fluxes,  in 
turn,  drive  the  evolution  of  the  mixed  layer,  the  inversion, 
and  thus  ducting  conditions  of  th®  lower  atmosphere. 


39 


The  052340Z  February  1980  sounding  was  used  to  '•es: 
model  sensitivity  tc  varying  subsidence  rates.  The  hindcast 
subsidence  was  -.  0030  m/s  (Table  II}  which  resulted  in  a 
predicted  395  meter  thick  elevat  ed  duct  with  the  top  at  740 
meters  elevation  and  the  base  at  345  Tatars.  Subsidence  was 
varied  from  this  baseline  value  and  deviations  in  the  pre¬ 
dicted  duct  height  arid  -hickness  were  recorded.  This  par¬ 
ticular  sounding  was  chosen  because  the  associated  initial 
and  verifying  ducts  were  elevated,  which  allowed  the  varied 
subsidence  tc  affect  the  position  of  the  predicted  duct  in 
both  the  upward  and  downward  directions. 

Deviations  in  the  3uct  predictions  resulting  from  subsi¬ 
dence  specification  errors  of  +/-  .0030  m/s  are  245  m  at  the 
top,  26  9  m  at  the  base,  and  25  m  in  thickness  (Table  VII). 
Errors  resulting  from  assuming  persistence  in  duct  heights 
(Table  III)  are  appr cximatel y  equal  to  the  errors  resulting 
from  a  subsidence  specification  error  of  +/-  .0030  m/s. 
From  this  single  case  comparison,  it  clear  that  the  model 
is  a  r;  im prove  men t  over  persisted  ce  only,  when  subsidence  can 
be.  specified  to  within  ♦ /-  .0030  m/s . 

The  range  of  subsidence  calculated  during  the  short 
period  from  6  to  22  February  1980  in  the  data  area  was 


-.  0030  to  -.0075  m/s  or  -.00525  +/-  .00225  m/s.  Assuring 
this  to  be  a  normal  range  or  variability  and  the  median 
value,  -.00525  m/s,  tc  be  representative  of  the  climat oiogi- 
cal  subsidence,  then  the  "dim  urological'1  value  appears  to 
be  less  than  the  +/-.0030  m/s  marginal  utility  criterion. 
This  conclusion,  based  or.  a  limited  data  set,  seems  to  indi¬ 
cate  that  a  climatological  subsidence  value  *ould  provide  a 
duct  height  forecast  which,  on  the  average,  has  some  skill 
ever  persistence  or  a  forecast  of  no  change. 

In  reali-y,  errors  will  not  be  confined  to  subsidence 
alone,  but  will  be  caused  by  errors  in  other  inputs  such  as 
wind  and  SST.  While  SST ,  wind,  and  subsidence  errors  may  be 
offsetting  at  times,  there  will  be  times  when  they  will  be 
additive.  When  all  the  errors  ere  additive,  the  expected 
duct  height  error  for  SST  plus  wind  measurement  errors  are 
163  m  at  the  top  and  154  m  at  the  bottom  (Table  V'i  .  Given 


•'•he 

erro 

rs  due  to  assuming  persist 

ence 

i:i  duct  heights  of 

364 

m  at 

the  •Lop  and  214  m  at  the 

base 

of  the  duct  (Table 

IV)  , 

the 

difference  of  approximately 

207 

n  at  the  -op  and  60 

m  at 

the 

base  indicate  the  maximum  i 

uct  h 

eight  error  permis- 

sible  due  tc  subsidence  specification  error  is,  in  the  mean, 
130  m.  From  Table  VI,  the  subsidence  specification  errer 


TABLE  VII 


Sens itivi 

t  y  cf 

Mo  del 

Duct  Pred 

ic tions 

to  Subside 

",Cc« 

SUBSIDENCE 

PREDICTS D 

DUCT  PARAMETERS 

(m  /s) 

TOP 

CHANGE  BASE 

CHANGE 

THICKNESS 

CHAi 

.0300 

-.0005 

’8*8 

2  80 

2  30 

645 

585 

300 

240 

37  5 

3  3  5 

-  20 

-10 

-.00  10 

910 

170 

535 

190 

37^ 

-20 

-.0015 

865 

1  25 

480 

135 

3  35 

-10 

-.00  20 

830 

90 

430 

85 

400 

5 

-.00  25 

785 

45 

380 

35 

405 

10 

-.00  30  •* 

740 

0 

34  5 

0 

34^ 

n 

-.0035 

713 

-27 

312 

-33 

401 

6 

-.0040 

675 

-  65 

250 

-85 

4  1  5 

20 

-.0045 

630 

-1  10 

225 

-120 

405 

10 

-.0050 

600 

-  1  40 

180 

-165 

4  20 

25 

-.0055 

570 

-1  70 

145 

-200 

425 

3  0 

-.0060 

535 

-2  05 

110 

-235 

425 

30 

EMS  3P.R0ES  : 

DEVIATION  ?  FOM  REFERENCE  EXPECTED  DEVIATIC 


SUBSIDENCE 

TOP 

BASE 

+/-  .0030 

245 

269 

♦/-  .0025 

202 

221 

+/-  .0020 

156 

178 

+/-  .0015 

1  18 

128 

+/-  .0010 

79 

85 

♦/-  .0005 

37 

34 

THICKNESS  (Ws  X  24  hours) 


25 

259 

22 

216 

23 

173 

10 

130 

15 

86 

8 
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*  reference  values 
**  heights  in  meters 
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TACTICAL  APPLICATIONS 


The  first  consideration  in  tactical  application  or  the 
model  is  the  recognition  of  its  capabilities  and  limita¬ 
tions.  It  can  be  used  to  estimate  when  low  clouds  or  fog 
can  be  expected  to  form.  It  can  be  used  to  estimate  the 
expected  position  of  low  elevated  and  surface  basec  ducrs 
below  1200  m  over  an  18  -  24  hour  period.  However,  it  can¬ 
not  predict  information  on  upper  level  ducts. 

The  tactical  significance  of  knowing  up  to  30  hours  in 
advance  that  visibility  may  be  reduced  or  a  low  ceiling 
developed  will  depend  on  current  and  planned  operations. 
The  results  presented  in  Chapter  IV  B  indicate  that  the 
model  does  well  in  predicting  low  level  stratus/fog. 

When  fog  or  low  level  stratus  clouds  are  predicted,  the 
following  possible  effects  on  own  and  enemy  forces  in  each 
warfare  area  should  be  considered: 

(1)  maneuverability  restricted 

(2)  lost  or  reduced  visual  signalling  capability 

(3)  reduced  visual  target  d  etection/i  lent  if  ication 

(4)  degraded  optical  equipment  performance 
including  IH  weapons  and  sensors 

(5)  hampered  flight  operations  for  carrier  based 
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aircraft  and  LAMPS 


Tactical  employment  of  Naval  assets  requires  a  knowledge 
of  the  state  of  the  environment  if  those  assets  are  to  be 
employed  effectively.  To  that  end,  IP.EPS  (Integrated 
Refractive  Effects  Prediction  System)  was  developed.  when 
radiosonde  data  are  input  to  IREPS,  a  refr activity  profile 
is  generated  which  is  utilized  to  define  the  location  of 
ducts  and  to  assess  radar  coverage,  ESM  and  communication 
ranges. 

Eecause  soundings  are  normally  taken  only  cnce  or  twice 
daily  from  aircraft  carriers  and  the  atmospheric  boundary 
layer  undergoes  change,  there  exists  a  need  to  be  able  to 
predict  changes  in  the  refract! vity  profile  between  sound¬ 
ings.  This  prediction  can  be  critically  important  to  battle 
croups  in  the  positioning  of  both  surface  and  air  assets, 
and  in  the  intelligent  management  of  the  EMCON  plan. 

The  duct  prediction  capability,  as  detailed  in  Chapter 
IV  A,  provides  a  significant  planning  tool.  Tactical  useage 
of  the  duct  predictions  stems  from  the  effects  of  ducts  on 
EM  propagation.  Prediction  of  a  surface  based  duct  indi¬ 
cates  probable  extended  ranges  for  transmitter-receiver 
antenna  pairs  in  the  duct  and  possible  holes  or  gaps  in 
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coverage  just  above  the  duct  foe  antennas  in,  near,  or  below 


the  duct. 

The  predicted  M  profile  may  be  used  as  an  environmental 
data  set  tc  run  IREPS  and  coverage  diagrams  generated  for 
specific  emitters.  This  procedure  was  used  with  a  hypothet¬ 
ical  surface  search  radar  under  observed  refractive  condi¬ 
tions  and  model  predicted  refractive  conditions.  IREPS 
version  1.7  (unclassified)  was  used  no  generate  the  coverage 
diagrams.  Radar  parameters  for  the  hypothetical  radar  were: 
Antenna  height  110  ft.,  antenna  type  sin(x)/x,  vertical  beam 
width  10°,  elevation  angle  0°,  free  space  range  55  nautical 
miles,  and  frequency  5  Ghz. 

Clearly  shown  in  Figs.  14  and  15  is  the  model's  value  in 
predicting  the  occurance  of  low  level  ducts.  The  radar 


HEIGHT  (ft) 


C  40K 


Figure  15.  Surface  Search  Radar  Coverage  0240  9  Feb  80 
Based  on  Actual  Sounding  (1)  and  Model 
Prediction  (r)  . 

(2)  radiosonde  data  are  input  to  IREPS  and  coverage 
diagrams  are  generated  to  make  current  assessments 

(3)  radiosonde  data  and  surface  observations  are  input 
tc  the  UPS  model  and  a  30  hour  prediction  is  made. 

(4)  *he  predicted  M  profile  is  input  to  IREPS  and 
predicted  coverage  diagrams  are  generated. 

(5)  operational  plans,  incorporating  forecast  ducting 
effects,  are  proposed. 

(6)  tactical  decisions,  using  the  environment  ^o  the 


advantage  of  the  force,  are  made 


In  a  scenario  similar  to  that  observed  in  the  Arabian 
Sea  7  to  9  February  1980  when  an  elevated  due*  became  sur¬ 
face  based,  many  tactical  considerations  are  involved.  Sup¬ 
pose  the  battle  group  wanted  to  avoid  detection  by  an  enemy 
force  in  the  region.  Forewarned  that  surface  based  ducting 
was  expected  to  occur  in  12  to  18  hours  and  that  (JHF  radio 
cransmissions,  normally  useful  between  ships  only  up  to  25 
*o  30  nautical  miles,  could  be  detected  at  ranges  in  excess 
cf  200  nautical  miles  due  to  ducting,  the  battle  group  com¬ 
mander  would  be  able  to  impose  a  more  restrictive  EMC0N  con¬ 
dition  in  adequate  time  to  prevent  hostile  intercept  of  the 
group's  transmissions.  Alternatively,  if  the  battle  group 
commander  expected  tc  be  opposed  by  active  surface  search 
radar,  he  could:  initiate  a  preemptive  attack  before  his 
presence  was  exposed;  steam  out  of  range;  or,  relax  EMCON  at 
the  appropriate  time. 

For  a  battle  group  in  a  multithreat  environment,  early 
threat  detection  is  critical  to  defeat  the  threat  at  minimum 
cost  to  the  battle  group's  offensive  capability.  The  12  to 
18  hour  advance  notice  of  surface  based  ducting  provides  the 
tactical  planner  the  •‘■ool  necessary  to  use  his  limited  sur¬ 
veillance  assets  with  the  efficiency  needed  for  sustained 
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operations  and  the  e ffective nas s  require!  for  early  threat 
detection.  Surface  ships  could  be  positioned  in  such  a  way 
as  to  take  maximum  advantage  of  the  extended  surface  search 
radar  ranges  afforded  by  the  surface  based  duct.  Airborne 
surveillance  assets,  freed  from  task  of  surface  surveil¬ 
lance,  could  concentrate  their  efforts  toward  detecting  and 
tracking  air  targets.  The  net  result  is  increased  force 
ef f ectiveress.  Depending  on  the  battle  group's  speed  and 
the  picket  station  assigned,  the  s  irface  ship  may  require 
several  hours  to  reach  the  assigned  station,  thus  making  the 
leal  time  provided  by  the  prediction  all  the  more  important. 

The  t ;  mine  of  war  at  sea  strikes  to  coincide  with  pre¬ 
dicted  occurrence  of  a  surface  based  duct  could  give  consid¬ 
erable  advante.ge  to  the  strike  aircraft  by  enabling  them  to 
use  predicted  holes  just  above  the  duct  to  penetrate  enemy 
surveillance  and  to  use  the  duct  for  targeting  and  jamming. 
In  such  a  case,  a  12  to  18  hour  prediction  provides  time  to 
clan  mission  tactics,  arm  aircraft,  and  conduct  mission 
briefs. 

When  sonobuoy  patterns  are  laid  and  monitored  by  LAMPS, 
one  constraint  on  placement  of  the  bouys  is  LAMPS'  on  sta¬ 


tion  endurance,  when  a  surface  based  duct  occurs,  soncbcuvs 


car:  be  monitored  at  greater  distances.  Prediction  i 
advance  of  surface  based  ducting  provides  the  lead  time  nec 
essary  tc  plan  and  lay  bouy  patterns  at  distances  from  th 
battle  group  which  take  advantage  of  the  extended  monitor!?. 
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•  £Q££i2SIp.tIS  AND  RECOMMENDATIONS 

The  Tactical  Environmental  Support  System  (TESS)  ,  cur¬ 
rently  under  development,  will  incorporate  various  satellite 
derived  products,  and  sinple  station  assessment  systems  such 
as  I3EPS  and  the  NPS  mixed  layer  model.  This  will  place 
state  of  the  art  environmental  sensing  and  predictive  capa¬ 
bility  at  the  disposal  cf  those  who  need  the  information: 
the  operational  commanders. 

Based  cn  evaluations  with  a  limited  data  set,  the  NPS 
model  has  been  demonstrated  to  perform  well  both  in  pre¬ 
dictions  of  lew  level  cloud  formation  and  in  predictions  of 
low  elevated  and  surface,  based  ducts.  Further  work  is 


required  to  determine 

t  '.vi  cn  o.i  a  1  •  s 

range  of 

applicab ilit 

both  geographically  and 

seasor.a  lly. 

It  is 

believed  tha 

when  subsidence  can  be  specified  with  ar.  accuracy  of 
+  /-.0015  n/s,  the  model  can  be  used  to  accurately  forecast 
low  level  ducts.  Two  additional  recommendations  are 
detailed  in  the  following  paragraphs. 

In  tactical  situations  where  the  outcome  may  depend  on 
accuracy  of  duct  predictions,  a  method  is  needed  to  verify 
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the  prediction.  A  radiosonde  launch  may  be  impractical  due 
to  EMCON  conditions.  An  acoustic  sounder  can  be  used  to 
measure  the  height  of  the  inversion  which,  as  shown  in  Fig. 
3,  corresponds  ‘o  the  top  of  the  duct.  Davidson  ,  et  al 
(1982),  while  conducting  marine  atmospheric  boundary  layer 
research,  utilized  the  economical  acoustic  sounder  on  R.  V. 
ACANIA  to  monitor  inversion  height.  Using  a  combination  of 
daily  radiosondes,  NPS  model  duct  predictions,  and  periodic 
acoustic  soundings  to  track  the  inversion  height  would  pro¬ 
vide  improved  capability  to  use  the  environment.  It  is 
therefore  recommended  that  similar  acoustic  sounders  be 
installed  on  aircraft  carriers. 

Determining  appropriate  temperature  and  humidity  pro¬ 
files  during  initialization  of  the  model  is  a  possible 
source  of  error.  There  are  two  schools  of  thought  which 
differ  on  the  appropriate  heiaht  for  the  inversion  initiali¬ 
zation.  Sounding  data  typically  show  the  inversion  as  a 
layer  100  to  200  meters  thick,  while  the  model  requires  a 
single  height  to  be  specified.  One  school  suggests  that  the 
bottom  of  the  inversion  layer  ie  best  while  ■‘■he  second  sug¬ 


gests  that  the  midpoint  of  the  layer  is  most  reasonable 


The  work  m  this  thesis  used  the  former  method.  The 
model  predictions  (Table  III)  show  a  bias  ir  that  nearly  all 
of  the  predicted  duct  heights  are  lower  than  these  observed. 
This  saems  to  indicate  that  inversion  heights  should  be 
digitized  at  the  midpoint  of  the  lay«r  instead  of  at  the 
bottom.  It  is  recommended  that  further  work  be  dene  in  this 
area  to  establish  a  mere  objective  method  of  digitization. 

The  NPS  atmospheric  boundary  layer  modal  could  provide 
the  fleet  a  predictive  capability  and  -actical  planning 
*ool.  Consideration  should  ba  given  to  introducing  it  into 
the  fleet  for  operational  evaluation.  The  need  to  under¬ 
stand  the  model  assump+ions  in  order  to  assess  model  appli¬ 
cability  in  a  given  environmental  situation,  the  need  for 
accurate  subsidence  estimates,  and  the  requirement  for  sur¬ 
face  wind  predictions  precluie  model  installation  except 
where  qualified  personnel  are  assigned.  In  this  regard,  the 
model  should  be  installed  at  locations  (aircraft  carriers) 
where  there  are  Geophysics  Officers  who  car.  interpret  and 
use  the  model  as  a  tool  along  with  other  meteorological 
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RADIOS  CNDE  TIMES  AND  LOCATIONS 


TABLE  VIII 
Radiosonde  Data  Sat. 


SOUNDING 

SHIP 

DTS 

LATITUDE 

LONGITUDE 

1 

MI  MITZ 

052340Z 

23 ° 1 2  ’  N 

Of.  10 09'  E 

2 

NIMITZ 

070020Z 

22  0  38  1  N 

061036*  E 

3 

CORAL  SEA 

070030Z 

17040*  N 

96  5°5  1 '  E 

4 

HI  MITZ 

0723  33Z 

2  3  0  3  7  *  H 

061000*  E 

5 

CORAL  SEA 

082240Z 

23°00  *N 

06 1 o  1 1  *  E 

6 

HI MITZ 

0823  33Z 

23  °1 3  *  N 

06 1 °00  *  E 

7 

NIMITZ 

1  100  14Z 

23° 1 7 1 N 

060059*  E 

8 

HI  MITZ 

1  12346Z 

2  3  °  1 0  *  N 

06 1 ®0  1 '  E 

Q 

HI MITZ 

132307Z 

23°03»N 

064°30*  E 

10 

NIMITZ 

1500  13Z 

2  3  o  1  1  «  n 

06 1 °07 ’  E 

11 

HI  MITZ 

2103  00Z 

24  0 0  1  *  N 

053032’ E 

12 

CORAL  SEA 

2120  35Z 

23°06  *N 

061°12'  E 

13 

NIMITZ 

2  12244Z 

2301  3*  N 

05  0058*  E 

14 

CORAL  SEA 

2220  15Z 

23°  3  3*  K 

060OQ9' E 

*  date-time  groups  February  1933 
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